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The formation of a finite density polaronic state is analyzed in the context of the Ifolstein model 
using the Dynamical Mean-Field Theory. The spinless and spinful fermion cases are compared to 
disentangle the polaron crossover from the bipolaron formation. The exact solution of Dynamical 
Mean-Field Theory is compared with weak-coupling perturbation theory, non-crossing (Migdal), 
and vertex correction approximations. We show that polaron formation is not associated to a 
metal-insulator transition, which is instead due to bipolaron formation. 
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Recent experiments strongly suggest that the electron- 
phonon (e-ph) interaction is relevant in many materi- 
als, ranging from high Tc superconducting cuprates I"], 
to colossal magnetoresistance manganites |2J, from the 
fullerenes to magnesium diboride 0. While the spe- 
cific role of the e-ph interaction is certainly very differ- 
ent from compound to compound, the properties of these 
materials are hardly explained in terms of standard ap- 
proaches to the e-ph problem, like the Migdal-Eliashberg 
theory, and pose a serious challenge to theories. 

More specifically, polaronic features have been ob- 
served in lightly doped cuprates 0, in the manganites 
up to some indication in the fullerenes A small 
polaron is a carrier so tightly interacting with the lattice 
that its effective mass is strongly enhanced therefore re- 
ducing its mobility!^. The single polaron problem has 
been extensively studied allowing to understand in detail 
the polaron physics [1,0]. Nevertheless, the experimen- 
tal findings of polaronic effects obviously deal with finite 
densities of carriers, and prompt for an analogous under- 
standing of the finite density polaron problem, where the 
polarons arc able to interact and to change drastically 
the phonon properties. In this regard, it is important to 
stress the distinction between polaronic and bipolaronic 
states. If repulsion between the electrons is neglected, 
two polarons (with opposite spin) tend in fact to bind, 
giving rise to a (local) pair, which is called bipolaron. 
The bipolarons in turn may undergo a superconducting 
transition of the Bose-Einstein type . If superconduc- 
tivity is not allowed, bipolarons give rise to an insulating 
state of localized pairs |llLll2| which may eventually con- 
dense in a charge-ordered state |0 . In some important 
compounds bipolaronic states are indeed unlikely formed. 
In cuprates and fullerenes, the strong electron correlation 
forbids bipolaron formation, while in the manganites, the 
double-exchange mechanism favors ferromagnetic states 
at low temperature, in which no bipolaron can be formed. 



In order to disentangle the polaron effects from bipo- 
laron formation, we compare the spinless fermion case, in 
which bipolarons are forbidden by Pauli principle, with 
the spinful case, which has been extensively studied in the 
recent past [3 E EE [13, El • To be more explicit, we 
show that the polaron crossover at finite density is not by 
itself a metal insulator transition (MIT), and insulating 
behavior can only be associated to localized bipolarons. 
This effect shows up in the non-vanishing of both the 
quasiparticle renormalization factor, and of the renor- 
malized phonon frequency in the spinless case. In the 
spinful case the quasiparticle weight vanishes at some def- 
inite coupling while the renormalized phonon frequency 
remains finite. 

We consider the Holstein molecular crystal model, in 
which tight-binding electrons interact with local modes 
of constant frequency. The Hamiltonian is 

H = -t c\ ^Cj,a+H.c.-g ni{ai+a\)-^uJo a\ai, 

(1) 

where Ci^^ (cJo-) a.i (a|) are destruction (creation) 
operators for fermions and for phonons of frequency cjqj 
Ui the electron density, t is the hopping amplitude, g is 
an e-ph coupling. The density is fixed ton = 0.5(n=l) 
in the spinless (spinful) system, which corresponds to the 
particle-hole symmetric half-filled case. In analogy with 
the studies of the Mott transition in the Hubbard model, 
in which antiferromagnetism is neglected fl9i|. we restrict 
ourselves to the state with no charge order [l^. Despite 
in the particle-hole symmetric case the ground state is 
always ordered, the homogeneous solution may become 
representative of the groundstate whenever charge or- 
dering is spoiled by some frustration effect like a next- 
nearest-neighbor hopping. Moreover, this study allows 
to characterize how strong e-ph interaction may lead to 
the destruction of the metallic state (just like electron 



correlation leads to the Mott state). The adiabatic ra- 
tio 7 = LOa/t has been shown to be an important pa- 
rameter for the single polaron formation. In the adia- 
batic regime 7 < 1, the polaron crossover occurs when 
A = /luqI ~ 1, while in the antiadiabatic regime 7 > 1, 
the condition is instead [g/u>oY — 1 [1)13 ■ 

We solve the model using the Dynamical Mean-Field 
Theory (DMFT), a non-perturbative approach which be- 
comes exact in the limit of infinite dimensions We 
notice that in such a limit the spinless fermion case does 
not coincide with the infinite correlation limit. Nonethe- 
less, this system represents an instructive playground 
where polaronic effects can be observed without too 
many competing phases. In DMFT, the lattice model 
is mapped onto an impurity problem subject to a self- 
consistency condition, which contains the information 
about the lattice. In our case the impurity model is 

k,(7 k,(T 

- gia + a^)J2Pj<r+u^oa^a, (2) 

(7 

where the phonons live only on the impurity (/) site, Ek 
and Vfc are the energy levels and the hybridization pa- 
rameters of the conduction bath. For the z-coordination 
Bethe lattice of half-bandwidth D = 2t^ — 1 the self- 
consistency equation in the z — > 00 limit is 

k 

In this work we use Exact Diagonalization (ED) to solve 
the impurity model Q 20] . This method requires to re- 
strict the sum in Eqs. (|2I3|I to a finite small number of 
levels Ns — I. The discretized model can then be solved 
at r = using the Lanczos method. The convergence 
is exponential in Ns, and a few levels are enough to give 
converged results. Results are obtained taking typically 
Ns = 10, having checked that no significant change oc- 
curs for larger Ns jl^. Error bars can be evaluated by 
linearly extrapolating (in 1/Ns) to A^s ^ cx), and are of 
the size of symbols in the figures. We compare the ED 
solution of DMFT, with two approximate schemes: a self- 
consistent Non Crossing Approximation (NCA) [ill Il5| 
(diagrams (a) in Fig. ^ and a self-consistent approxima- 
tion including the first Vertex Correction beyond NCA 
(VCA) [ill (diagrams (a) and (b) in Fig. [TJ. Notice 
that in the standard Migdal-Eliashberg approximation 
the phonon spectrum is not self-consistently evaluated, 
but it is taken "from experiments" , while in our NCA 
and VCA the phonon self-energy is self-consistently eval- 
uated (on the imaginary frequency axis) through an iter- 
ative scheme, just like the electron self-energy. Actually, 
to allow for a safer truncation on the frequency axis, the 
quantity which is iteratively determined is not the full 
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FIG. 1: The NCA [(a)] and the VCA [(a) + (b)] diagrams 
for electron (upper graphs) and phonon (lower graphs) self- 
energies. Wavy lines are phonon propagators while straight 
lines are electron propagators. Second order perturbation the- 
ory is obtained replacing internal lines with free propagators. 

self-energy, but the difference between this quantity and 
the result from second order perturbation theory. Zero 
temperature results are obtained by lowering the tem- 
perature until the physically relevant quantities converge 
(consistently increasing the number of Matsubara freque- 
cies Nmax = 6/(2T7r)). We checked the convergence to 
T = by monitoring the quasiparticle spectral weight 
Z defined as Z"! = (1 - (S(iw„=i) - EiiuJn=o)/2TTT) 
where a;„ — {2n — 1)t:T. The vanishing of Z is used 
within DMFT to characterize the Mott MIT in the re- 
pulsive jl9j, and the pairing transition in the attractive 
Hubbard models [2^ • A vanishing Z has been found also 
in the spinful Holstein model [13. The convergence to 
T — turns out to depend on both 7 (as detailed below, 
we study 7 = 0.1 and 7 = 1 as representative of adiabatic 
and nonadiabatic regimes, respectively) and A. Within 
VCA T/t = 3 X 10~^ is sufficient to get results repre- 
sentative of the ground state for weak/intermediate cou- 
pling. In the spinful case for 7 = 0.1 T — 10"^ is instead 
necessary since the polaron crossover is approached for 
smaller coupling (see below). Within NCA it is possible 
to span the strong coupling regime, making an extrap- 
olation to T = necessary in the adiabatic regime for 
A > 1 (A > 0.5) in the spinless (spinful) case. In the non- 
adiabatic regime the extrapolation is required for A > 1.6 
(A > 1.2) in the spinless (spinful) case. 

We first discuss the spinless case. Exact DMFT results 
for Z are shown in Fig. El The logarithmic scale on the y- 
axis evidences that in both cases Z, even if exponentially 
reduced, does never vanish by increasing A, indicating 
that no MIT is taking place. It is important to observe 
that Z increases when the truncation error is reduced in- 
creasing Ns ■ The comparison with NCA and VCA shows 
non trivial tendencies. Both approximations are accu- 
rate at weak coupling, but VCA remains closer to ED 
for relatively large coupling even in the adiabatic regime, 
where the Migdal approximation would be expected to 
hold. Quite surprisingly even in the non-adiabatic regime 
VCA improves NCA only at weak coupling. At strong 
coupling NCA predicts a polaronic crossover, in qualita- 
tive agreement with exact results, even if, for 7 = 1, the 
crossover coupling is strongly underestimated by NCA. 
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FIG. 2: Spinless fermions. The quasiparticle weight Z within 
second order perturbation theory (2PT), NCA, VGA, and ED 
for 7 = 0.1 (upper panel) and 7 = 1 (lower panel) as a func- 
tion of A. The arrow in the upper panel marks the MIT for 
7 = P3 . The inset shows ED and NGA in a larger range 
of A for 7 = 1 to make the polaron crossover visible. 



On the other hand, VGA drastically diverges from ED, 
and gives unphysical negative mass renormalization at 
strong coupling both for 7 = 0.1 and 1, in agreement 
with the divergence of vertex corrections predicted in Ref. 
[m . The coupling at which the approximate methods de- 
viate from exact results decreases with increasing 7, and 
it is unrelated with the polaron crossover coupling which 
is instead larger for 7 = 1, as shown in the inset of Fig. 
1 , and in agreement with the single polaron case 0, 1^ ■ 

An important difference between the finite-density sit- 
uation and the single polaron problem is that many elec- 
trons are able to renormalize the phonon properties. The 
renormalized phonon frequency can be obtained from 
the phonon propagator I?(ia;„) = J^^ dre'^'^"'^ {T(ao{T) + 
alir)){ao + al)) , as {fl/ujo)^ = -2/ujQD{iujr, ^ 0)-^. As 
it is shown in Fig. |21 fl/uiQ never vanishes as a function 
of the coupling, but rather exponentially decreases. The 
way NCA and VGA results for fl compare with ED is 
analogous to the results for Z. In the adiabatic regime 
NGA strongly overestimates ft/ujo around the polaron 
crossover. In the nonadiabatic regime NGA again pre- 
dicts a phonon softening which does not actually occur 
at such small values of the coupling. In the same regime 
VGA improves NGA result only at weak coupling and 
gives a phonon hardening at strong coupling. It is worth 
to note that for both Z and the coupling at which 



FIG. 3: Spinless fermions. The renormalized phonon fre- 
quency fl/iOQ. Notations are the same of Fig. |21 



NGA and VGA deviate from ED is lower in nonadiabatic 
than in adiabatic regime. This "non uniform" behavior 
is not present for a single polaron 9], and represents a 
first peculiarity of the finite-density case. 

Now we compare our findings for spinless fermions with 
the spinful fermion case. In the adiabatic regime a MIT 
has been found around A ~ Xmit — 0.68 close to 
the value (A = 0.664) at which the density of state at 
Fermi energy vanishes in the adiabatic limit (7 = 0) |l6| . 
This transition has a precursor in the phonon softening 
but, contrary to the spinless case, the spin degrees 
of freedom prevent electron coherent hopping at strong 
coupling through a Kondo-like mechanism 1231. A MIT 
has been claimed to occur also within NGA |ll| . 

In Fig. 21 we show results for 7 — 0.1. The logarithmic 
plot clearly shows that Z vanishes faster than exponen- 
tial, signaling a MIT at A ~ 0.76, in agreement with 
Ref. 12]. Gontrary to the spinless case, here Z decreases 
by increasing Ng. On the other hand, Q/loq decreases 
with A as in the spinless case, suggesting that, even if 
strongly softened, the phonon mode never becomes com- 
pletely soft. This behavior is more evident for larger 
phonon frequency 7 = 1, as reported in Fig. Also 
in this case Z vanishes much faster than exponential at 
A = 1.44, where the MIT takes place, while the phonon 
renormalization is much less effective and leads to quite 
a large value of fl/uQ even at the MIT point (see inset). 

The comparison with approximate schemes qualita- 
tively resembles the spinless case. NGA underestimates 
Z for A below the MIT, but it gives a finite weight also 
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above the exact MIT point. Our extrapolation to T = 
from finite temperatures does not allow us to definitively 
rule out the existence of a MIT within NCA even if, in 
contrast with a claim in Ref . (ll| , Z seems to decrease ex- 
ponentially with A within this approximation |24l | . Again 
VGA gives better results than NCA even in the adiabatic 
regime up to intermediate coupling A ~ 0.55, but it be- 
comes completely unreliable at strong coupling. 

We have studied the formation of a polaronic state in 
the Holstein model at half- filling within DMFT. For spin- 
less fermions, a continuous crossover leads to a polaronic 
state by increasing the coupling constant. Despite the 
electron effective mass becomes exponentially large in the 
strong-coupling regime, the ground state is always metal- 
lic. The crossover is more abrupt in the adiabatic case. 
In the spinful case, the polarons can bind to form bipo- 
larons, leading to a real MIT. The phonon renormaliza- 
tion is much stronger in the adiabatic regime than in the 
nonadiabatic case, but the phonons do not become com- 
pletely soft at the MIT. Approximate treatments (NCA 
and VCA) strongly deviate from exact DMFT above a 
coupling which diminishes with increasing u^^jt. At weak 
coupling VCA correctly reproduces the qualitative trends 
of exact results, and improves on NCA. 
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FIG. 4: Spinful fermions. Z and fi/ojo (inset) for 7 = 0.1 as 
a function of A. The arrow marks the MIT for 7 = p^. 



[1] 
[2] 

[3] 
[4] 



A. Lanzara et al, Nature 412, 510 (2001). 

G. Zhao, K. Conder, H. Keller, and K.A. Miiller, Nature 

381, 676 (1996). 

O. Gunnarsson, Rev. Mod. Phys. 69, 575 (1997). 
J. Kortus, et al, Phys. Rev. Lett. 86, 4656 (2001). 



N 



0.1 r 1 



0.01 



0.001 




FIG. 5: Spinful fermions. Z and Q/uio (inset) for 7 = 1.0 as 
a function of the coupling. 



[9: 

[10 

[11 
[12: 

[is: 

[14 

[is: 

[16 
[17 

[is: 

[19 
[20: 
[21 



[22: 
[23: 

[24: 



p. Calvani, et al, Phys. Rev. B 53, 2756 (1996) 
M. Ricco, et al, Phys. Rev. B 67, 024519 (2003). 
Here we will only be concerned with small polarons aris- 
ing from Holstein-like interactions. For the definition of 
other kind of polarons see, e.g., J. T. Devreese, in Ency- 
clopedia of Applied Physics (VCH, Weinheim (1996)). 
M. Capone, W. Stephan, and M. Grilli, Phys. Rev. B 56, 
4484 (1997); S. Ciuchi, F. de Pasquale, S. Fratini, and D. 
Feinberg, ibidem 56, 4494 (1997). 

M. Capone, C. Grimaldi, and S. Ciuchi, Europhys. Lett. 
42, 523 (1998). 

A. Alexandrov and J. Ranninger, Phys. Rev. B 23, 1796 

(1981) 

J. P. Hague and N. d'Ambrumenil, |cond-mat/0106355| 
D. Meyer, A.C. Hewson, and R. Bulla, Phys. Rev. Lett. 
89, 196401 (2002). 

S. Ciuchi and F. de Pasquale, Phys. Rev. B 59, 5431 
(1999). 

J. K. Freericks, M. Jarrell, and D. J. Scalapino, Phys. 
Rev. B 48, 6302 (1993); J. K. Freericks, ibidem 48, 3881 

(1993) . 

J. K. Freericks, Phys. Rev. B 50, 403 (1994) 

A. J. Millis, R. Mueller, and B. I. Shraiman, Phys. Rev. 

B 54, 5389 (1996). 

P. Benedetti and R. Zeyher, Phys. Rev. B 58, 14320 
(1998). 

A. Deppeler and A. J. Millis, Phys. Rev. B 65, 100301 
(2002). 

A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, 
Rev. Mod. Phys. 68 13 (1996). 

M. Caffarel and W. Krauth, Phys. Rev. Lett. 72, 1545 

(1994) . 

We also truncate the phonon Hilbert space introducing a 
cut-off in the number of phonons Nph, and always check- 
ing for convergency. The truncation is more delicate in 
the adiabatic regime, where a large number of phonons 
can be excited. When necessary, we used Nph ~ 50. 
M. Capone, C. Castellani, and M. Grilli, Phys. Rev. Lett. 
88, 126403 (2002). 

G. Moeller, Q. Si, G. Kothar, M. Rozenberg, and D. S. 

Fisher, Phys. Rev. Lett. 74, 2082 (1995). 

The exponential decay of Z within NCA shown in Fig. 0| 



hardly identified using a linear scale as in Ref. [Tl| . 



